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Future sattelite (MAP and Planck), balloon and ground based experiments will 
determine the basic cosmological parameters within a few percent. We examine 
here the effect of this on constraining the SUSY parameter space for supergravity 
i?-parity conserving models (with tan f3 < 25) for the cases of fCDM and ACDM 
cosmological models. For the vCDM (ACDM) models, the gluino mass is restricted 
by m g ~ 720(540)GeV. In both cases, the cosmological constraints are sensitive to 

non-universal SUSY soft breaking producing a lower bound mj ^ AOOGeV in some 
regions and for the i/CDM model, gaps in the allowed rrig range for other regions. 
For gluino (neutralino) masses greater than 450(65)GeV, mo is constrained to be 
small making squark and slepton masses generally light and determined mostly by 



1 Introduction 

One of the significant phenomena that has arisen with the development of 
supersymmetry has been the deepening of the connection between particle 
physics and cosmology. While this connection has always been present (e.g. 
from the early nuclear synthesis calculations El) supersymmetry (SUSY) now 
offers models that reach up in energy to the GUT scale Mq rs 10 16 GeV (and 
perhaps further to the Planck scale) and backwards in time to the very early 
universe. In particular, models with i?-parity invariance automatically predict 
the existence of dark matter, and for a resonable range of SUSY parameters, in 
an amount comparable to what is observed astronomically. These models have 
been waiting to be tested, and there are now a large number of experiments 
that are currently on line or will be on line in the relatively near future that 
will determine whether supersymmetry is a valid theory of nature. 

While up to now the flow of information has been mostly from particle 
physics to cosmology, the new sattelite experimentsa MAP and Planck as well 
as about 25 balloon and ground based experiments can reverse this and can 
give constraints on what might be expected at accelerators. These astronomi- 
cal experiments will measure the cosmological parameters with great accuracy. 
One uses the parameter Qi = pij p c where pi is the density of matter of type 
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1 and p c = 3H 2 /8itGn (H =Hubblc constant, Gjv =Newton's constant) to 
measure the mean amount of matter of type i in the universe. Current mea- 
surements give H = h (100 km s _1 Mpc^ 1 ) with 0.5 t h * 0.75 (and hence 
p c = 1.88 x 10 - 29 h 2 gm / cm 3 ) . For a number of cosmological models one finds 
for cold dark matter (CDM) the range 

0.1 < n C DAih 2 < 0.4 (1) 

and this range of CDM is the one assumed in many particle physics calcu- 
lations. However, MAP and Planck measurements will greatly restrict this 
window, and impose strong constraints on the SUSY parameter space, and in 
this way restrict the SUSY mass spectrum. Thus astronomical measurements 
will impinge upon what is expected to be seen at accelerators. 

2 Supergravity Models 

In order to analyse dark matter predictions we need a well-defined SUSY 
model. We will use here supergravity grand unified models with gravity medi- 
ated SUSY breaking (at a scale Z Mq) in a hidden sector, and with gravity 
being the messenger of SUSY breaking to the physical s£C±orl3. This gives 
rise to the following types of soft breaking masses at M^M: too (scalar soft 
breaking mass), TO4/2 (gaugino soft breaking mass), Aq (cubic soft breaking pa- 
rameter) , Bo (quadratic soft breaking parameter) , as well as non-universal soft 
breaking termsu. In addition, a Higgs mixing parameter, po, for the two Higgs 
multiplets Hi and H 2 {W {2) = ^H^) can form naturally with p,Q « 

In the following we will assume that the gaugino masses are universal at 
Mq (as would be true to a good approximation for a simple GUT group). 
Over much of the parameter space p 2 ^> Af§ which leads to the "scaling 
relations" for the charginos (x^,i = 1,2) and neutralinos = 1 — 4) 

2to x o £ m x ± £ m x o S (| - \)m~ g . 

We assume here that the first two generations of squarks and sleptons are 
degenerate at Mq (to suppress flavor changing neutral currents) with common 
mass too and parametrize the Higgs and third generation soft breaking masses 
at Mq as: 

m 2 Hl = m 2 {l + 5^); m 2 H , 2 = m 2 (l + S 2 ) (2) 

m\ L = ml(l + S 3 ); m 2 UR = m%{l + Si); m\ R = m%{\ + S 5 ) (3) 

m\ R = mo(l + «5a); mf L = to 2 (1 + 8 7 ) (4) 

Here qL = (uL,rfi) is the L squark doublet, ur the R up-squark singlet, 
etc. We limit our parameters to the domain mo,TOg < 1 TeV; \A t /mo\ < 
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7; tan/3 < 25; \5i\ < 1 where A t is the i-quark A-parameter at the electroweak 
scale. These represent a choice of "naturalness" conditions. Note that for 
tan /3 'Z 25, £5, 5q, 67 make only small contributions and can be neglected. 
(They would become important for larger tan/3.) 

The RGE0 then give for /i 2 at scale Mz the result a: 
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Here D = 1 - (rn 4 /200 sin/3) 2 , A fl ^ A f - 0.613mg, S = TjrFm 2 (Y = 
hypercharge, m 2 = masses at Mq) and Cg is given in Ibahez et alQ Do vanishes 
at the t- quark Landau pole and hence for m t = 175 GeV is generally small: 
Dq < 0.23. (Ar is the residue at the Landau pole.) Note that the choice 
62 < 0, Si > and S3 + S4 > will increase the size of fi 2 , while 82 > 0, Si < 
and £3 + £4 < will decrease the size of /i 2 . In the following we will see that n 2 
plays a key role in dark matter predictions, and the effects of these particular 
sign possibilities of Si will allow a qualitative understanding of the phenomena. 



3 Direct Detection of Dark Matter 

We review briefly in this section the direct detection of dark matter particles 
in the Milky Way incident on a terrestial detector. The density of such matter 
is estimated at pdm — 0.3 GeV/ cm 3 with an impinging velocity of vdm — 
300 km/s. Calculation of detector event rates proceeds in two steps. One 
first calculates the expected Xi relic density of f2 x o h 2 , left after annihilation in 
the early universe i. The size of ^ x «h 2 varies as one moves across the SUSY 
parameter space. In particular, for m x o < 65 GeV (or by the scaling relations, 
nig t 450 GeV) the annihilation cross section ^dominated by the s-channel 
Z and h poles and requires sensitive treatment E^fEd. For m x o ^ 65 GeV, the 
t-channel squark and slepton poles become dominant. These two regimes will 
show up below in the detector event rates. 

One first restricts the SUSY parameter space so that fl x o h 2 falls within the 
allowed window of the cosmological model under consideration. In addition, 
one restricts the parameter space so that accelerator bounds are obeyed. One 
then calculates the scattering of incident Milky Way \i by quarks in a nuclear 
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target □ within the above restricted parameter space, giving rise to detector 
event rates R measured in events /kg d. 

4 Cosmological Models 

Current astronomical measurements have sufficient uncertainty that they al- 
low for a large variety of cosmological models. However, the Planck and MAP 
sattelites will be able to determine the basic cosmological parameters, H, A, f2 
etc. to (1-10)% accuracy by measurements of the deviations AT from the CMB 
temperature To = (2.728 ± 0.002) °K. Different cosmological models predict 
different angular correlations of the CMB deviations allowing for experimen- 
tal determination of the cosmological model. We consider here two possible 
cosmological models that might result from such measurements, and examine 
what consequences these might have on the sparticle spectra, and hence on 
accelerator searches for supersymmetry. 

(i) vCDM Model. If neutrinos have masses of O(eV) they could furnish the 
hot dark matter (HDM) of the universe, f2„. In addition there is baryonic dark 
matter (B) and CDM (the neutralinos x?) such that the total £1 = 1. As an 
example we assume that the sattelite measurements determine central values 
to be fl^ = 0.2, £Ib = 0.05, SIcdm = 0.75 and a Hubble constant of h = 0.62. 
Using estimates of the errors with which these quantities can be measured by 
the Planck sattelite E3o we find 

VcDMh 2 = 0.288 ±0.013 (6) 

This shows the remarkable accuracy future determinations of cosmological pa- 
rameters are capable of when compared with the current knowledge of Eq. 

©• 

(ii) ACDM Model. As a second model we assume the existence of a cos- 
mological constant A along with baryonic and CDM, and assume that Planck 
has determined the central values of these parameters to be = 0.05, 
&CDM = 0.40, f^A = 0.55._and h — 0.62. Using the estimates of the errors in 
the Planck measurements E3 one finds here that 

VL CDM h 2 = 0.154 ±0.017 (7) 

While the above models are only two possible examples, they represent bounds 
on what might actually exist. The important point is that the future balloon 
and sattelite measurements will be able to determine these parameters well, 
and the above models represent limits within current knowledge. 
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vCDM(lo), 5,=0=5 2 
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Figure 1: Maximum and minimum event rates for a xenon CDM detector vSjijluino mass 
for fCDM model (1 std window) for universal soft breaking, fi > 0. From Rcftj. 

5 vCDM Model 

We consider first universal soft breaking where by Eqs. (g)-(Q), the 8i = 0. Fig. 
1 plots the maximim and minimum event rates for Xe detector as a function 
of the gluino mass [nig = (7 — 8)m x o by the scaling relations] for a 1 std band 
of Eq. (§). One may compare this with Fig. 2 of i where only the broad 
band of Eq. (|lj) is imposed. One sees a reduction in event rates in the region 
nig * 450 GeV (i.e. m x o < 65 GeV), the Z and h pole dominated region in 
the early universe annihilation, and a corresponding increase for the higher 
masses. More significant is the appearance of forbidden regions, i.e. gaps, in 
the allowed values of wig in the region m g = 500 GeV and m g = 600 GeV. In 
examining the effects of non-universal soft breaking, we note the correlation 
that when /i 2 is decreased (which by Eq. occurs for 82 = 1 = —81) the event 
rate R is increased and when /i 2 is increased (82 = — 1 = — Si) R is decreased. 
Fig. 2 shows these effects for the latter case, and also that the forbidden 
regions below m g = 500 GeV are significantly widened. For the alternate 
possibility (82 = 1 = —81), the event rates increase. The maximum rates for 
this case are at the sensitivity of the current Nal detectors E3, and detectors 
with a sensitivity of R Z 10~ 3 events /kg d would cover the entire parameter 
space. The gaps of the previous case have now disappeared. However, m g is 
now bounded from below i.e. m g Z 420 GeV. In all of these cases there is 
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vCDM(la), S 1 =l=-8 2 
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Figure 2: Same as Fig.l for non-universal soft breaking with <5i = 1 = S 2 . From RefQ. 

also an upper bound on rrig of rrig t 720 GeV. Thus cosmological constraints 
bound the allowed gluino mass range, and are sensitive to the non-universal 
soft breaking. 

The large i-quark mass causes large L — R mixing in the stop (mass) 2 
matrix, making the light eigenvalue, to- , small. In general mj is governed 
by the size of too and in the light ncutralino domain too can be large, i.e. 
mjj can rise to w 600 GeV . However, for heavy neutralinos (m^ ~ 65 GeV, 
rrig Z 450 GeV) the early neutralino annihilation is governed by the t-channel 
sfermion poles, and in order to get sufficient annihilation for Eq. (^) to hold, 
Too must be small i.e. Too t 200 GeV. One finds in this domain that mj is 
generally an increasing function of m g with mass ranging from 250 GeV to 
500 GeV. 



6 ACDM Model 

Here fl x o h 2 obeys Eq. (Q) which even more sharply restricts the SUSY param- 
eter space. Fig. 3 shows the maximum and minimum expected event rates 
for a Xe detector for different values of Si. While one has the same general 
behavior as in the j/CDM model, there are interesting differences that distin- 
guish the two. Thus for this case for the 1 std window of Eq. (]?]), one finds 
a much tighter upper bound of rrig < 520 GeV. There is also a lower bound 
rrig ^ AOOGeV for <5 2 = 1 = —Si (dashed curve), sharply restricting the gluino 
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ACDM(la) composite 
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Figure 3: Maximum and minimum event rates for a xenon CDM detector vs. gluino for 
ACDM model (1 std window), /i > 0, for universal soft breaking (solid curve) and non- 
universal soft breaking <5i = 1 = — 82 (dotted curve), <5i = —1 = — 82 (dashed curve). (The 
discrete high mass points are for 81 = — 1 = —82 where scaling is badly broken.) From 
RefR 



mass range for this case. Also, there are no gaps in the allowed gluino mass 
range. 

The squarks and slcptons are again constrained to be relatively light in the 
domain rrig Z 420 GeV since here Eq. (Q) requires mo ^ 100 GeV. Thus mi 
ranges from (200 — 400) GeV in this domain and for the first two generations 
of squarks, rriq ranges from 400 GeV to 500 GeV and approximately scales 
with nig. The lightest slepton, en can have a mass as low as (85 — 90) GeV in 
this region, at the edge of detectability of LEP 190. 

This work was supported in part by National Science Foundation Grants 
PHY-9722090 and PHY-96020274. 
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